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We show that a proper consideration of the contribution of Trugman loops leads to a fairly low 
effective mass for a hole moving in a square lattice Ising antiferromagnet, if the bare hopping and 
the exchange energy scales are comparable. This contradicts the general view that because of the 
absence of spin fluctuations, this effective mass must be extremely large. Moreover, in the presence 
of a transverse magnetic field, we show that the effective hopping integrals acquire an unusual 
dependence on the magnetic field, through Aharonov-Bohm interference, in addition to significant 
retardation effects. The effect of the Aharonov-Bohm interference on the cyclotron frequency (for 
small magnetic fields) and the Hofstadter butterfiy (for large magnetic fields) is analyzed. 
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I. INTRODUCTION 

Ever since the discovery of high-temperature super- 
conductivity in cuprates,^ understanding the motion of 
a hole in an otherwise half-filled Cu02 layer has been a 
major challenge in condensed matter physics.- If. as most 
customary, one models this system with a one-band Hub- 
bard model,' then in the limit of a large on-site repul- 
sion U the half-filled case maps onto an antiferromagnetic 
(AFM) Heisenberg Hamiltonian with spin-exchange cou- 
pling J = At'^/U, where t is the nearest-neighbor (NN) 
hopping. ' In this limit, then, the problem reduces to un- 
derstanding the motion of a hole in a two-dimensional 
(2D) AFM background.^"^ 

A major reason for the difficulty in dealing with this 
question is that despite being long-range ordered at 
T = 0, the undoped AFM background has a very compli- 
cated wavefunction due to spin fiuctuations, which make 
it very unlike the semi-classical Neel state. A simple but 
more realistic description that could be used for analyt- 
ical purposes is missing; as a result, progress has been 
made primarily through numerical simulations. '"'" 

Here, we present an essentially exact numerical solu- 
tion, and an approximate but quite accurate analytical 
approach, for the much simpler issue of a hole moving in 
an Ising AFM background,'^" whose undoped wavefunc- 
tion is the Neel state. The solutions have a variational 
interpretation, and are shown to be accurate for a wide 
range of parameters, of up to t/ J ^ 3. 

The generally accepted view, which probably explains 
why this problem has not been solved so far (to the best 
of our knowledge), is that a hole cannot really move in 
an Ising AFM, since as it hops away from its initial loca- 
tion it rc-shufflcs the spins at the sites it visits, creating 
a string of "defects" (misaligned spins). In dimensions 
higher than one, the energy cost of this string increases 
roughly linearly with its length, and as a result the hole 
is forced to stay in the vicinity of its original position, i.e. 



its effective mass is infinite.'' In this view, spin fluctua- 
tions which effectively remove (or "heal" ) pieces of this 
string of defects are needed to free the hole and allow it 
to acquire a finite effective mass. Of course, such spin 
fiuctuations are absent in an Ising AFM. 

That the hole is not truly localized even in an Ising 
AFM has been known ever since Trugman pointed out'' 
that by going twice around a closed loop, the hole can 
actually acquire a finite effective mass. We will return 
to this in more detail below, but the main idea is that 
while the first circuit along the closed loop creates the 
usual chain of defects, the second reshuffle of the spins 
during the second circuit removes all these defects, but 
also ends with the hole at a different location than the 
original site. By repeating such processes, the hole can 
therefore move anywhere on its original sublattice. 

Even though Trugman started from a Neel background 
and identified this mechanism for generating a finite 
quasiparticle mass in the absence of spin fiuctuations, he 
included spin-fiuctuations in his calculation by consid- 
ering the full Heisenberg AFM Hamiltonian when deter- 
mine the effective mass of the hole. ' The issue of the hole 
mass in a purely Ising AFM thus remained unanswered. 

Here we carry out this calculation, and show that the 
hole is fairly mobile if i ^ J. Moreover, if a transverse 
magnetic field is turned on, due to Aharonov-Bohm inter- 
ference of the Peierls phases associated with these closed 
loops, the effective hoppings acquire a magnetic field de- 
pendence over and above the usual Peierls phases, with 
interesting consequences. Our results reveal that even 
this seemingly simple problem is actually very interest- 
ing and leads to rather non-trivial results. 

The paper is organized as follows. In Sec. II, we specify 
the model and our notation. In Sec. HI, we describe the 
analytical calculation in the smallest variational subspace 
where the hole acquires mass, both with and without 
magnetic field. We then explain the generalization for the 
numerical calculation. Section IV contains our results, 
and a summary and final conclusions appear in Sec. V. 



II. MODEL 



Consider a spiii-i Ising AFM on a square lattice 






jY^ {cr^,zCTj,, ~ I) , (1) 



where cFz is the Pauh matrix, and j = J/4. The ground 
state is the classical Neel state, with all spins on sublat- 
tice A up, and all spins on sublattice B down. 



iGS)=n4tndiio)> 



(2) 
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where Ci^a annihilates an electron at site i, and for con- 
venience we shifted the energy so that "HafmIGS) = 0. 

We would like to study the dynamics of a single hole 
in this system, within the approximation that only NN 
hopping of electrons is possible, however no-double occu- 
pancy is allowed. As a result, at each site we have either 
the hole, or a spin. Moreover, the spin can only be in its 
proper orientation (consistent with the sublattice its site 
belongs to) or flipped (i.e., a frozen magnon-like "defect" 
is created at this site). 

In the following wc will only keep track of the location 
of the hole and of the flipped spins. To simplify the 
notation, wc introduce a "defect" creation operator 



Hi e A 
Hie B 



(3) 



where a are the raising/lowering Pauli matrices, and 
hole creation operators 



Ci.-|-, if i G A 
Cj4, Hi e B 



(4) 



for a hole at the site i. With this notation, for example 
hid AGS) means that the hole is at site i and the spin at 
site j y^ i is flipped; all other sites have their spins in the 
proper Neel configuration. 

Consider now the motion of the hole, with the no- 
double occupancy condition enforced. If the hole is at 
site i, the only possibility is for an electron from one of 
its four NN sites j to hop into i, thus moving the hole to 
site j. If the spin of the electron at j was properly ori- 
ented, when it moves to i it has the wrong orientation; in 
other words, a "defect" is created at i when the hole hops 
from i -^ j. On the other hand, if there was a "defect" 
at j to begin with, when the electron moves to i it will 
be properly oriented, and therefore the "defect" at j is 
removed as the hole hops from i ^- j. 

Thus, the Hamiltonian that describes the dynamics of 
the hole in the 2D Ising AFM is 



n = p 






h\h, 



7t 



H.c. 



P + Uafm- (5) 



The projector P enforces the no-double occupancy con- 
straint, as well as the condition that each site has either 
the hole, or a spin: h\hi+d\di + did\ ~ 1. In the presence 
of a uniform transverse magnetic field B, the hopping in- 
tegrals tij include the proper Peierls factors (see below). 

This Hamiltonian is similar to the Edwards model, ^^ 
however, unlike that model it enforces (i) the fact that 
there can be no magnons at the site where the hole is, 
i.e. states like /i]d||GS) are forbidden; (ii) the fact that 
there can be at most one spin-flip per site, i.e. states 
like /i^(d])"|GS),n > 2 are forbidden. Finally, (iii) un- 
like the Edwards model, where the energy of the defects 
was assumed to be described by ^^idldi, here wc use 
"Hafm to calculate the true cost for creating spin fiips - 
two neighboring defects cost less energy than two farther- 
apart ones, because they only disrupt seven AFM bonds, 
not eight. Such corrections are less important than is 
imposing (i) and (ii), but keeping track of the proper 
exchange energies is simple enough and we do so. 

Note that a description of a Heisenberg AFM would 
require the addition of terms oc 2j^,^ j\{didj + djdp, 
describing the XY spin-exchange interaction leading to 
spin fluctuations. However, since the Neel ground-state 
|GS) is not a good approximation for the undopcd ground 
state of this model, it is hard to quantify the meaning 
of adding these terms. Spin fluctuations could also be 
introduced with a small magnetic field parallel to the 
X-axis, leading to a term oc X]i('^I + '^Oi ^■s considered 
in the Edwards model.' ' As shown there, such terms do 
significantly lower the effective mass of the hole, but we 
ignore them here. 



III. FORMALISM 

A. Hole propagation when B = 

We begin by considering the case where no transverse 
magnetic field is applied. In this case, the unit cell con- 
tains two sites (one from sublattice A, one from sublat- 
tice B) and the corresponding magnetic Brillouin zone is 
a square rotated by 7r/4 (see below). For each k in this 
Brillouin zone, we define the plane wave 



^E 






(6) 
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and want to calculate the Green's function 



G(k, uj) = (GS|/ikG(w)/ij^|GS) , (7) 



where 



G(c.) 



1 



cj + it] — H 



(8) 



is the resolvent associated with this Hamiltonian, rj > 
is infinitcsimally small, and we set h = 1. N denotes 



<»' □ D n n D □ □ 
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FIG. 1. (Color online) Shortest sequences of repeated hops 
around closed loops leading to effective 2nd NN (top panel) 
and 3rd NN (bottom panel) hopping. The hole is marked 
by a blue square, and the defects (spin flips) it creates or 
annihilates as it moves are marked by red circles. 



the number of unit cells, or the number of sites in each 
sublattice, and is taken to infinity. 

Besides producing the single-hole spectrum from its 
poles, the spectral function A(k, oj) = — ilmG(k, w) is 
the quantity that would be measured by spin-polarized 
ARPES, assuming that a (photo)electron with spin-up is 
removed from the system. Of course, there is a second set 
of wavefunctions consisting of plane waves involving the 
sites on sublattice B; their associated Green's function 
would be related to ARPES measured when a spin-down 
electron was ejected from the system. Obviously, the 
spectral weights are identical for the two cases. 

As discussed, as the hole hops away from its initial site, 
it creates a string of defects as it reshuffles the spins on 
its path. These defects increase the energy of the state 
roughly linearly with the length of the string, because 
of the broken AFM bonds. In the absence of spin fluc- 
tuations, there are only two ways to remove this costly 
string of misaligned spins. One is for the hole to retrace 
its path, removing all the defects - in this case it ends 
up at the original site it started from.^'^ Such terms can 
only renormalize the overall energy, but do not gener- 
ate dispersion. The second possibility is for the hole to 
go repeatedly around closed loops: on the first circuit a 
string of defects is created, while on the second circuit 
the defects are removed. These loops generate effective 
second and third NN hopping terms, as counted in the 
original square lattice. 

In Fig. 1 we show examples of the shortest sequences 
of events that result in generating such 2nd (a) and 3rd 
(b) NN hoppings. ' One can easily convince oneself that 
longer closed loops can only contribute to one of these two 
effective hoppings, since the final location of the hole is 
within two hops of the original one. Effective NN hopping 
is impossible, because of spin conservation: if the hole is 
on the other sublattice, there must be an odd number of 
magnons around. 

In perturbational terms, for t <^ J — Aj, the sequence 
of hoppings depicted in Fig. 1 (a) results in 2nd NN hop- 
ping i2 = t''/(6^ • 10^ • 12/), since the initial and final 
states, which are both in the no-defect manifold, are con- 
nected through 6 hopping processes, and the five inter- 
mediary states have energies higher by 6J, lOj, 12j, lOj, 6J 



respectively. Similarly, Fig. 1(b) gives a contribution to 
3rd NN hopping tg = t^°/{6^ ■ 10^ • 14^ • 18^/), as the 
longer strings of defects are more costly. 

Note that clockwise hopping of the hole around the 
same loop as shown in Fig. 1(b) would generate a higher- 
order contribution to the 2nd NN hopping tj ~ t^ be- 
cause some of the intermediary states have slightly differ- 
ent energies. This allows us to estimate the range where 
perturbation theory is valid. If we ask that ^2/^2 < 0.1, 
so that this correction from 5-defect processes to the 3- 
defect result is small, we find t < 8.4J = 2.1J. In other 
words, because of the relatively long and thus energeti- 
cally costly loops of flipped spins involved in generating 
the effective hoppings, perturbation theory extends up to 
rather large values t < 2J. For comparison, in cuprates 
it is believed that t ^ J/3, so they would be just outside 
the perturbational range (for these values, ^2 ~ 0.4^2). 
Of course, the effects of spin fluctuations cannot be ig- 
nored in the cuprates. 

Clearly, then, for t<2J the dominant contribution to 
the hole's dynamics comes from the processes of Fig. 1(a). 
We therefore begin by calculating G(k, ui) in a variational 
approach, by only allowing such short loops (with up to 
3 defects on 3 corners on a square plaquette) to be gener- 
ated. Longer loops are more expensive and the probabil- 
ity to generate them should be lower. Of course, they can 
be included in a numerical calculation, but this will lead 
to only quantitative, not qualitative changes (see below) . 

The method we use is similar to that used in Refs. 13 
and 14 for the ID, and in Ref. ..■> for the 2D Edwards 
model. Within this variational approximation, we gener- 
ate the equations of motion for the propagator by using 
repeatedly the Dyson identity 



G{lj) ^ Goiio) + G{Lu)VGoiio) , 



(9) 



where Go{uj) is the resolvent for Ho = "Hafm, while V is 
the first term in Eq. (5), describing the hopping of the 
hole. We then find 



G(k, w) - Go(w - 4j) [1 - i ^ Fi (k, w, 5) 



(10) 



where Gq{uj) = l/{Ljj+iT]), 4j is the cost of having only the 
hole in the system (first state in Fig. 1(a)), and d points 
to any of the four NN sites, i.e. S £ {(±1,0), (0,±1)} for 
a lattice constant a = I. The new propagators are 



Fi(k,w, 






■JkRi 



{GS\huG{Lo)dlhl^s\GS) , 



ieA 



(11) 

and are related to the amplitude of probability to propa- 
gate between a no-defect state and a 1-defect state such 
as the second state shown in Fig. 1(a), with a defect at 
the original site. In the following, for simplicity we will 
denote i^i(k, w,<5) as Fi{S). 

An equation of motion for Fi can now be generated. 
The hopping V can lead to three different outcomes: ei- 
ther the hole hops back to site i, removing the defect; or 



it hops by another S, creating a hnear string of the type 

'^l'^+sH+2s\^^)'^ O'^ ^t hops by one of the two d' _L S, 
leading to a state such as the 3rd state shown in Fig. 1(a). 



Only the first and last outcomes are allowed within our 
variational calculation, leading to 



F,{6,) = -tGoi^ - lOj) 



G{k,L,)+ Y, ^2(^1,^2) 



where 






N 



(12) 



(13) 



describes generalized Green's functions associated with 2-defect states. Its equation of motion, within our variational 
space, is 



where 



^2(^1,^2) - -tGo(c^ - 14 j) [F,{Si) + ^3(^1,^2,-^1)] 



ieA ^jV 



(14) 



(15) 



describes states with a string of three consecutive defects. Only states such as shown in the middle panel of Fig. 1(a) 
arc in our variational space, hence the unique F3 term in the equation of -F2. Finally, the equation for Fg, is now 
connected to two possible 2-defect states. The hole could either hop back, linking to -^2(^1,^2), or it could complete 
the closed loop by hopping to site i and removing the defect that was there, resulting in 



„jkRi 



J2^^{GS\h^G{c.)hldl,^dl,^^JGS)=e-^^('^+'^^F,{^S,,-6,) 



(16) 



Eq. (16) follows after a translation by —(^1 + ^2), which maintains the site i on the same original sublattice A, since 
82 -L ^1. Of course, hopping also links F^ to two generalized propagators F^ with 4 defects, however those are ignored 
in this 3-defect variational space. Thus, we find: 



F3iSi,S2,-di) = -tGo{u-mj)[F2{Si,52) + e-'''<''+'-^F2{-S2,-Si] 

In other words we have a closed system of linear equations for G and the various allowed Fi^ F2, F3 functions. 
This system can be solved analytically. The final result is 



(17) 



G(k,w) 



1 



cj + i?/ — eo('^) + 4f2('^) cosfcj; coaky + 4t3(aj) [cos(2fca;) + cos(2fcj,)] 

I 



(18) 



where 



eo(w)=4j + 4t2Go(w)-4t3, 
2t2Go(w)/3(a;) 



LJ = UJ — lOj 



t2iuj) = — 



l-4/32(cj) ' 



tsiuj) = /3{u;)t2{uj) . 
The functions that appear in these definitions are 

e-Go{<:j)Go{i^Mu;) 



Pico) 



1-jHuj) 



(19) 

(20) 
(21) 

(22) 



2t2Go(J)) 



and 



-f{Lj)^t^Go{d;)GoiLo^l6j), 



w = w - 14j - t'^Goiuj - 16J) . 



(23) 
(24) 

(25) 



The Green's function of Eq. (18) describes a free par- 
ticle on the A sublattice, with an onsite energy eg, and 
hoppings t2 and t3 that correspond to NN and 2nd NN 



fi 



FIG. 2. (Color online) An effective fs process generated with 
only 3-defect strings. 



for sites on this sublattice. Interestingly, all these effec- 
tive quantities are functions of w, in other words retar- 
dation effects are explicitly taken into consideration in 
this variational calculation. For t ^ j we find that in- 
deed eo « 4j- 4iV6j, while fs = <V(2f600/). This is 
precisely twice the value we estimated for Fig. f (a); the 
factor of 2 is due to contributions from both clockwise 
and anticlockwise loops. 

The appearance of the t^ term at this level of the cal- 
culation is surprising at first, since configurations with 
5 defects like that appearing in Fig. 1(b) are not al- 
lowed in this variational space. The answer for how to 
achieve t^ hopping with only 3-defect loops is shown in 
Fig. 2: essentially, instead of removing the last defect 
at the 6th hopping, another 3-dcfect process is initiated. 
It is straightforward to check that in the perturbational 
limit, the ratio between the expected values for this t^ 
and the ^2 of Fig. 1(a) is of i''/(7200j^). This is precisely 
the ground-state value of /3(aj = 4j) in this asymptotic 
limit, which further validates this statement. It is also 
now clear that this strategy of starting a new loop when 
only a single defect is left can be repeated an arbitrary 
number of times, and an infinite sequence of higher order 
contributions to ^2 and t^ can be generated this way. This 
explains the 1 — 4/3^(a;) denominators in Eqs. (20), (21), 
showing that our variational calculation sums contribu- 
tions from these 3-defect processes to all orders. 

As argued previously, the contribution of 5-defect loops 
is negligible if f < 2 J = 8j. Interestingly, for such values 
2/3 ~ i'^/(3600j^) varies from to just above 1 (however, 
this ignores the change in (3{uj) because the ground-state 
energy moves away from cj = 4j, as t/J increases). In 
other words, the denominators in Eqs. (20), (21) could 
become very small, implying that significant hopping, 
and thus light effective masses, may be possible in such 
systems even for rather small t/J ratios. 

As a final comment, the fact that k-dependence ap- 
pears only in the equation for F^ confirms that the effec- 
tive hopping terms are due to completely circling around 
the closed loops. If we truncated the variational space to 
include only 2-defect states, which would mean setting 
-Fj — > in Eq. (14). the result would be just a renormal- 
ization of the on-site energy but no dispersion. 



B. Hole propagation when B ^ 

Since the dispersion of the hole is due to Trugman 
loops, it is interesting to consider what happens when 
a transverse magnetic field B is applied. Because of the 
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FIG. 3. (Color online) Magnetic unit cell in real space when 
(j> — 27rp/(2n), where p and n are mutually prime integers. 
The unit cell is marked by the red dashed line, and contains 
2n sites, labeled by squares, on the A sublattice. These are 
indexed by their corresponding a — 0,...,2n — 1. The value 
of a also appears in the Peierls phases for hopping in the 
direction of the arrow. 



large number of elementary hopping processes involved in 
generating effective longer ranged hoppings, one would 
expect the phases associated with these effective hop- 
pings to be different from the normal Peierls phase. Put 
it another way, one would expect Aharonov-Bohm-like 
interference between clockwise and counterclockwise con- 
tributions, over and above the usual Peierls phases, and 
therefore a spectrum that is not simply that of a bare 
particle with hoppings t2,i3 placed in a magnetic field. 

Let $ = Ba^ be the magnetic flux through the ele- 
mentary square plaquette. Because the dressed hole only 
lives on one sublattice (sublattice A in our calculation), 
the actual relevant flux is 2Ba^, since the unit cell for 
the sublattice is doubled in size. 

From elementary considerations of the Hofstadter 
butterfly,^'' we know that we need to consider the cases 
when the ratio between the magnetic flux and the ele- 
mentary flux $0 = h/e is 2_Ba^/$o — p/^i where p and 
n are mutually prime integers. In this case, we expect the 
spectrum to split into n magnetic subbands. Of course, 
then, the Brillouin zone is folded down n times. We are 
therefore interested in the spectrum of the dressed hole 
when 



$0 



P_ 
2n 



(26) 



(for simplicity of notation, we also introduce (f> = 
27r$/$o = np/n). 

We use the Landau gauge A(r) = Bxey, in which 
case only the hopping integrals in the j/-direction acquire 
Peierls phases. The magnetic unit cell will contain 2n 
sublattice A sites, and it is drawn, together with the 
Peierls phases, in Fig. 3. Each site in the unit cell is la- 
beled by the index a = 0, . . . ,2n~l which establishes its 
location in the unit cell, as well as the phase associated 
with hopping off this site in the positive y-direction. 

The Hamiltonian is identical to that of Eq. (1), with 
the proper phases included in the hoppings. Of course, 
we assume that the magnetic field is not sufficiently large 
to overcome the AFM coupling and favor a FM ground 
state. Also, we ignore the Zceman term - since the spins 



are reshuffled as the hole hops around, the total z-axis 
spin Sz = X^i '^i,!^ remains constant and therefore the 
Zeeman contribution is an overall constant. 

We now need to introduce 2n distinct plane waves, 
associated with each site of the unit cell 



^Q,k 






(27) 



where A^ collects all sites on sublattice A with the same 
value of a. The procedure to generate the equations of 
motion for the propagator 



G„.„(k,c^) = (GS|V,kG(L.)/ilklGS> 



(28) 



is just as before, however now one has to take into consid- 
eration the Peierls phases. These can either be included 
in the definition of the new i^i,_F2,-F3, functions, or can 
be explicitly pulled out. In either case, the overall struc- 
ture of the system of linear equations is very similar to 
that in the B = Q case, and can be solved similarly. This 
allows us to reduce it to a single equation involving only 
G functions, which reads: 



G„'a(k,a;) = Go(w-eo) {(J^'a - <2 [2cos(fcy- (a- i)0)e''''="Ga',a-i(k,cj) + 2cos(fc.y - (a + i)^)e '''==G„',a+i(k, w)] 

-h [e2"=-Ga.,„_2(k,w)-F2cos(2fcy-2a</))Ga.„(k,w)-^e-2»'=-Ga'.a+2(k,w)]} . (29) 

Here we used the short-hand notations: 

4i2Go((i)(l-2/32(,.)) 



eo = eo(^,'^) =4j- 



1 - 4/?2(a;) -I- 4/34 (cj) sm^{2(t)) 



(30) 



h = t2{4>,u:) = -- 



2i2Go(w)/3(w) (1 - I3'iuj)) cos f + /32(w) cos ^ 



50 



1 -4/32(0;)+ 4/34(0;) sin^ (20) 



and 



is = t-i{4>,uj) = - 



2t^Goi^)l3^{uj)cos{2,p) 
1 -4/32(0;)+ 4/34(0;) sin^ (20) 



(31) 



(32) 



Remarkably, Eq. (29) is identical to the equation of 
motion one would obtain for the propagator associated 
with a Hamiltonian on the sublattice A only, with an on- 
site energy eg and NN and 2nd NN hopping (as defined 
on the sublattice) given by ^2,^3, plus the proper Peierls 
phases for the applied magnetic field. The terms in the 
parenthesis multiplied by ^2 come from NN hopping by 
ztSx ± Sy, with the proper Peierls phases accounted for 
in the phases of the cosine functions; similarly, the terms 
in the parenthesis multiplied by t^ come from next-NN 
hopping by i:2dx, ±2Sy. Only the ±2Sy hopping accumu- 
lates a Peierls phase, in this case, hence the cosine term 
appears only for that term. 

It follows that £0,^2 and is are the effective parame- 
ters describing the new quasiparticle (the dressed hole in 
the Neel AFM background plus magnetic field). They 
depend strongly on the magnetic field, besides the w- 
dependence due to retardation effects. As expected, 
they have the correct values when — > 0, as given by 
Eqs. (19)-(21). 

The cos -^ and cos -^ terms in ^2 , respectively cos 20 



r 



term in ia, are due to Aharonov-Bohm interference be- 
tween clockwise and counterclockwise loops. This can be 
checked directly in the asymptotic limit i ^ J by count- 
ing the accumulated Peierls phases for the processes from 
Figs. 1(a) and 2. Because of these, we expect the result- 
ing Hofstadter butterfly to have an unusual periodicity. 
That is, we still expect the appearance of n bands when 
2<I>/$o = p/n, due to interference between the Peierls 
phases. However, additional dependence of the hoppings 
on because of the Aharonov-Bohm interference will fur- 
ther modulate the overall bandwidth with a period 2tt in 
(the band structure is invariant if only ^2 changes its 
sign). In contrast, normally, (i.e. for constant £0,^2,^3) 
the spectrum is an even function of with period tt, 
because if the magnetic fiux through the unit cell (2$, 
in this case), increases by $0j the spectrum stays un- 
changed. Thus, we expect a quite different-looking Hofs- 
tadter butterfly for the dressed hole. In particular, since 
at least for small t/J we expect t^ <^ t2, the spectrum 
should become extremely narrow every time ^2 ^^^ be- 
cause of destructive Aharonov-Bohm interference. 



(a) 






(b) 




FIG. 4. (Color online) Examples of discarded configurations: 
(a) the central 6-defect configuration is only allowed to be 
linked back to the 5-defect configuration that generated it, 
shown to its right, or to the 7-defect configuration shown to 
its left. The other two "self-crossing" 5-defect configurations 
shown in (b) are not kept in the variational space. 
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C. Numerical calculation method 

The numerical calculation that we describe in this sec- 
tion has been carried out only for the case B = 0. This 
is a straightforward generalization of the 3-defect varia- 
tional calculation discussed above. The idea is to system- 
atically increase the size of the variational space, to see 
whether addition of configurations with longer loops of 
defects changes results significantly. If it does not, then 
wc know that the calculation is converged and therefore 
the results are essentially exact. 

We use an index TV to characterize the maximum num- 
ber of defects allowed within the variational subspacc. 
The equations of motion are generated as before, starting 
from just the hole and keeping all configurations linked 
to by hole hopping, up to those involving N defects. The 
exception is for longer chains with A^ > 4 defects, where 
we do not allow the chain to "self-cross" itself. For in- 
stance, for the 6-defect chain sketched in Fig. 4, the hole 
is not allowed to hop either left or down. Either of these 
processes would result in a configuration with the hole 
having two strings of defects attached to it. Trying to 
remove all these defects involves a very complicated se- 
quence of Trugman loops, and is therefore expected to 
have a very small contribution to ^2 and ts. The other 
role of such configurations is to renormalize energies, but 
as we will see in the next section, we can argue that their 
effect must be quite negligible. 

This is why for a configuration like that of Fig. 4, we 
allow the hole only to hop back (towards right) to the 
i^5 configuration from which it derived; or to hop up to 
generate a Fi configuration. In its turn, hopping to the 
left from this Fj configuration will close the loop and 
give a new contribution to the effective <2- New possible 
closed loops (and therefore additional contributions to 
^2, ^3) appear for A^ = 3, 5, 7, . . . 

Once all the allowed configurations have been gener- 
ated for a given N, the resulting (very sparse) linear 
system comprising their equations of motion is solved 
numerically using PARDISO.'^' We show results with up 
to iV = 7, in the following, because here there are al- 
ready over ten thousand allowed configurations, and, as 
we will see, this is sufficiently large for convergence to be 



FIG. 5. (Color online) (a) GS energy vs. t/J, as predicted 
by the analytic solution (circles) and numerically in the vari- 
ational spaces with A'' = 3,..., 7; (b) and (c) show the de- 
pendence on N of the GS energy for t/J = 3, respectively 5. 



achieved up to a quite high t/J ratio. 

Finally, note that in the numerical calculation, A^ = 3 
contains many more configurations than in the A^ = 3 
analytical calculation, because now all possible 2- and 
3-defect configurations (not just those like in Fig. 1(a)) 
are included. The results of the two methods, therefore, 
should not be expected to be identical. 



IV. RESULTS 

A. B — case 

We begin by looking at the ground-state (GS) energy 
as a function of t/J, using both the analytical and the 
numerical methods. Results are shown in Fig. 5(a) for 
t/J < 3. As expected, in the limit t/J -^ 0, all curves 
converge towards Eqs = J = 4j, which is the cost of 
placing the hole in the lattice (breaking four AFM ex- 
change bonds). As t increases and the hole acquires a 
finite mass, its energy is lowered. For small t/J, all val- 
ues are in excellent agreement, but for larger t/J they 
begin to fan out. At t/J = 3, convergence has already 
been achieved for A^ = 7, as shown in Fig. 5(b). For 
t/J~5,N = 7is not yet fully converged, although fur- 
ther corrections are not expected to be large. Note that 
we do not show results up to t/J = 5 in panel (a) be- 
cause the corresponding increase of the energy range to 
be displayed makes it even more difficult to distinguish 
the numerical curves from one another. 

The analytical calculation indeed gives a reasonably 
accurate value up to t/J ^ 1 ... 2. The fact that the 
biggest variation is between it and the A'^ = 3 numerical 
calculation, means that the 2- and 3-defect configura- 
tions ignored in the analytical approximation play a sig- 
nificant role in renormalizing the overall energy. This is 
not surprising, since these are some of the least expensive 
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FIG. 6. (Color online) Same as in Fig. 5, but for the GS 
quasiparticle weight instead of GS energy. 



configurations. However, qualitatively and even quanti- 
tatively it is clear that this analytical approximation is 
quite good up to fairly large t/J values. 
The quasiparticle weight 



Z{k) = \{^u\hl\GS)\' 



(33) 



defined as the overlap between the lowest eigenstate at a 
given momentum, 



H|(/.k> = i?(k)|0k> 



(34) 



and the non-interacting state /iJ^jGS), is shown at k = 
in Fig. 6. It remains quite considerable even for large 
t/J, showing that a significant part of the wavefunction 
consists of the hole alone, with no strings of defects. The 
effect of increasing the variational space is more visible 
here, as expected due to normalization: as the wavefunc- 
tion acquires extra components in a larger variational 
space, the weight of the hole-only part decreases, even if 
overall the energy of the state is not much changed. 

The hole band dispersion £'(k) and quasiparticle 
weight ^(k) are shown in Fig. 7, along high-symmetry 
cuts in the Brillouin zone of the original square lattice, 
for t/J = 3,iV = 7. The band folding due to the AFM 
order is clearly apparent. These results are very similar 
to those found for comparable parameters in the 2D Ed- 
wards model (see Fig. 10 of Ref. 15): the minimum is at 
k = (0,0), and k = (^, ^) is a saddle point (unlike in 
cuprates, where it is the actual minimum). The quasi- 
particle weight is fairly constant but with local variations 
that mimic the shape of the dispersion. This similarity 
is not surprising; even though the Edwards model allows 
more states (e.g. with more defects at a site), those are 
higher in energy and do not contribute much to the GS. 
The mechanism for generating an effective mass is iden- 
tical in both models, with only small quantitative differ- 
ences due to the energy assigned to the strings of defects 
in the two models. 

One difference between the two models is illustrated 
in Fig. 8, where we show the dispersion for TV = 3,5,7 
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FIG. 7. (Color online) Energy _B(k)/i (upper panel) and 
quasiparticle weight Z(k) (lower panel) along the high- 
symmetry directions of the square lattice Brillouin zone, for 
t/J = 3 and N — 7. The inset in the bottom panel shows 
the upper right quadrant of the full Brillouin zone; the AFM 
Brillouin zone is shaded. Lines are guides to the eye. 
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FIG. 8. (Color online) Dressed hole dispersion _B(k)/i along 
high-symmetry cuts in the Brillouin zone, for f/J = 3 and 
Af = 3, 5, 7. Symbols are numerical data (for N = 7, same 
as in Fig. 7) and lines are results to a fit as in Eq. (35), with 
parameters extracted for best fit on the (0,0) — (tt, 0) cut. 
Dashed lines are fits to extract the effective mass. 



and t/J = 3 along fewer cuts. The symbols are data 
from the numerical simulations, and the lines are fits to 
an effective dispersion 



- 4t2 cos fca; cos fcj, — 2^3 [cos(2fc2,) + cos(2fcj,)] , 



i?*(k)-i?o 

(35) 

along the (0, 0) -^ (0, tt) line. Using the same parame- 
ters along the (0, 0) -^ (tt, tt) line is clearly not a good fit. 
In contrast, for the Edwards model such fits worked well 
in the entire Brillouin zone. This suggests that the re- 
tardation effects may be somewhat stronger in this case: 
even though we know that only t2{yi) and ^3(0;) effective 
hoppings are generated, these may not be well approx- 
imated by constant values f 2 ■, ^3 if the lo dependence is 
considerable within the bandwidth of the dressed hole. 



This is less of an issue at smaller t/ J ^ where the band- 
width is significantly narrower. In fact, even for t/J = 3, 
the results for A^ = 3 shown in Fig. 8 clearly have a 
smaller bandwidth than for N = 5,7. This is not surpris- 
ing, since longer loops with additional contributions to 
effective hoppings are included in the latter cases. Even 
though the decrease in bandwidth is not that large, it 
is clear that the fit is better for the TV = 3 case. This 
may explain why the fit worked well for the 2D Edwards 
model, '^^ where the solution was restricted to TV = 3 and 
the bandwidths were narrower also due to more costly 
defects (equivalent to larger J) than in Fig. 8. 

Since we cannot extract meaningful t2i^3 values, we 
instead calculate the effective mass 



d^E(k) 



(J f£^ CJ Ky 



(36) 



IkKo 



with fits as shown by the dashed lines in Fig. 8. As 
expected, the mass is found to be isotropic. The results 
are shown in Fig. 9, on a logarithmic scale. The effective 
mass m* is in units of the band mass m, i.e. the bare 
particle mass if there was no AFM background. 

We show results up to t/J = 5, even though the val- 
ues are not fully converged for t/J > 3. At these larger 
values, longer loops than N = 7 would need to be in- 
cluded. From results such as in Fig. 5(c), however, we 
do not expect those further corrections to be very signif- 
icant. Also, they would decrease m* , as they would open 
additional channels for effective hopping, so the values 
shown in Fig. 9 can be taken as an upper bound for to*. 

In the limit t/J — >■ 0, the effective mass becomes very 
large. As discussed, from perturbation theory here we 
expect ^2 = t^/(21600j^). The corresponding m* /m = 
10.55( J/i)^ is shown as a dashed line in the inset, in good 
agreement with the values obtained from the analyti- 
cal and numerical calculations (symbols and full lines). 
The effective mass diverges as t/J —!• since the hole is 
bound to its original site in this case. As t/J increases, 
the effective mass decreases significantly, and it reaches 
m* /m ~ 40 for t/J = 3. Since this decrease is mirrored 
by the analytical iV = 3 calculation, it must be due to 
the summation over repeated 3-defect loops, like those 
shown in Fig. 2. Longer loops further lower the effective 
mass. The N = 5 loops give a significant contribution for 
t/J > 1, validating the expectations based on perturba- 
tion theory. The contribution of the N — 7 loops is still 
quite small, for these values. We expect to* to further 
decrease with increasing t/J — > oo. However, capturing 
that limit is very difficult if not outright impossible with 
this variational formulation, given the huge increase in 
possible configurations with increasing N. 

The main result, thus far, is that a hole in an Ising 
AFM is fairly mobile if t and J are comparable, even 
though there are no spin fluctuations in this model. 
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FIG. 9. (Color online) Effective mass vs. t/J. The inset 
shows the same data, plus the prediction of the weak-coupling 
perturbation theory (dashed line). 



B. B / case 

We now proceed to discuss the spectrum of the dressed 
hole in the presence of a transverse magnetic field B. We 
will use the A^ = 3 analytical Eq. (29), and solve it nu- 
merically for various values of ^ = np/n. We show results 
for t/J ~ 3, even though we know that here the analytical 
calculation is not sufficient for full convergence, simply 
because its S = bandwidth is sufficiently large to make 
it easier to see the effect of a finite B. We note that 
the numerical calculation in the variational spaces corre- 
sponding to various A^ can be carried out as well, however 
each previous unknown, corresponding to an allowed con- 
figuration of defects, now becomes a 2?! x 2n matrix of 
unknowns, corresponding to each of the 2n distinct sites 
in the magnetic Brillouin zone shown in Fig. 3. This huge 
increase in the size of the linear system to be solved, espe- 
cially for larger n values (smaller magnetic fields) makes 
the implementation of this scheme cumbersome, and un- 
necessary since we do not expect any qualitative changes. 

After solving the linear system in Eq. (29), we calculate 
the local density of states (LDOS) 



p{uj) = Im V'Gaa(k,a;), 



(37) 



where the sum is over the (highly folded) magnetic Bril- 
louin zone corresponding to the magnetic unit cell of 
Fig. 3. This LDOS is independent of the site a = 
0, ..., 2n — 1 used, hence the lack of a site index. 

We first consider the effect of the additional depen- 
dence on B due to the Aharonov-Bohm interference, on 
the Hofstadtcr butterfly expected when the magnetic flux 
through the unit cell is comparable to $0- Results for 
various simple ratios are shown in Fig. 10. In all cases, 
the upper panel shows the full results, whereas the lower 
panel shows the results when the Aharonov-Bohm in- 
terference is turned off, i.e. we use the (f> — values 
eo(0 = 0,a;), ^2(0 = 0,a;) and ^3(0 = 0,a;) in Eq. (29). 
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FIG. 10. 



(Color online) Local density of states p{i^) 



Loj J for various magnetic fluxes through the sublattice unit cell 



0, 



and 1, for i/J — 3. The upper panels show the full solution of Eq. (29), while the lower panels show the 



solution if the Aharonov-Bohm interference is turned off by setting (j!> = in Eqs. (30)-(32). The broadening is ry/J = 2.5 • 10 
except for the orange (light) curves in the insets, for which 77/ J — 2.5 • lO^**. 



For the 5 = case shown in (a), the results are of 
course identical. The density of states resembles that 
of a 2D square lattice with NN hopping. This is ex- 
pected, since the dominant hopping term is i2, which 
plays the role of NN hopping for the sublattice on which 
the dressed hole lives. The LDOS is slightly distorted due 
to (small) contributions from the tz term and the retar- 
dation effects. These small distortions and asymmetries 
are observed in all other panels. 

For a finite flux through the sublattice unit cell 
2$/$o = l/'T-, the band splits into n subbands, as ex- 
pected in standard Hofstadter butterfly phenomenology. 
The bands are not fully separated because we used a 
broadening 77/J = 2.5 • 10~^ which, while small, is still 
comparable with some of these features' bandwidth. De- 
spite this, the various subbands are easily identifiable. 

In Figs. 10 (b) and (d), the full result (top panel), while 
quite similar in aspect to the lower panel, shows a some- 
what narrower bandwidth. This is not surprising, since 
the Aharonov-Bohm interference terms like cos(3(?!)/2) 
etc. are responsible for a decrease in the value of the ef- 
fective hoppings, and therefore of the overall bandwidth. 
This also explains the apparent "collapse" of the butter- 
fly for (\) = 7r/3 and cj) = tt, shown in Figs. 10 (c) and (e), 
respectively. The very narrow peak seen in both cases 
shows the expected subband structure, if a much smaller 
77/J = 2.5 • 10~* values is used, as done in the inset: the 
former case shows the 3 subbands (these features are so 
narrow that even this much smaller ry can only partially 
resolve them), while the later case shows the one band. 
The extreme narrowing is due to the Aharonov-Bohm 
interference, which at these values of the flux results in 
a very small or vanishing t2{(j),i-u), see Eq. (31), and a 
bandwidth set by t3{(p,uj) <C ^2(0,^). 

As already noted, this additional B-dependent mod- 
ulation of the Hofstadter structure, coming from the 



Aharonov-Bohm interference, is also responsible for an 
increase in the periodicity of the butterfly. The lower 
panel shows that, as expected in models with constant 
hopping integrals, the butterfly is periodic if the flux 
through the unit cell increases by a flux quantum (here, 
2$ ^ 2$ -I- $0, or ^ (j) + tt). The full result clearly 
does not have this periodicity. As discussed previously, 
based on Eqs. (30)-(32) we expect the periodicity to be 
(j) ^ (j) + 2Tr for the full Aharonov-Bohm case; our simu- 
lations confirm this (not shown). 

We can imagine that this 27r periodicity survives if one 
allows longer loops in the calculation, but some care is 
needed. For example, loops like those shown in Fig. 1(b) 
would bring in factors of cos(70/2) in their contribution 
to ^2, and various other Aharonov-Bohm phases will be 
associated with longer loops. We calculated several of 
these and all are consistent with the 27r periodicity; so 
we believe this 27r period is correct to all orders, but do 
not have a proof. 

Finally, we discuss the role of the Aharonov-Bohm in- 
terference at very small magnetic fields. In this case, one 
expects the spectrum to separate in a sequence of Landau 
levels (LLs) separated by the cyclotron frequency. The 
appearance of LLs in the spectrum leads to quantum os- 
cillations in various transport measurements, such as de 
Haas-van Alphen oscillations. 

The LDOS at the lower edge, with (full lines) and 
without (dashed line) very small magnetic fields applied, 
is shown in Fig. 11(a). The B = LDOS shows the 
jump at the band edge, here smoothed out by a finite 
T] value, and the roughly constant 2D density of states 
above it. The slight monotonic increase is due to the use 
of a tight-binding dispersion, as opposed to an approxi- 
mate parabolic one.^'^ When the small field is turned on, 
we see that this LDOS splits into equally spaced narrow 
peaks, marking the LLs. The cyclotron frequency, de- 
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FIG. 11. (Color online) (a) LDOS for small magnetic fields 
(full lines), showing several LLs at the bottom of the spec- 
trum. The dashed lines shows the B = LDOS; (b) the lo- 
cation of the lowest LL peak, in the full calculation (squares) 
and with the Aharonov-Bohm interference turned off (circles). 
The dashed line is the analytical prediction of Eq. (38). These 
results are for t/J — 3 and 77/ J — lO^''. 



fined by the distance between consecutive LLs, increases 
roughly linearly with B, and so does the weight in each 
LLs, as expected because of its larger degeneracy. 

In Fig. 11(b), we plot the energy of the lowest LL 
against the flux through the unit cell. Squares show the 
full result, while circles are the result when the Aharonov- 
Bohm oscillations are turned off. The dashed line shows 
the value of the lowest LL for a particle of constant mass 



£=[^GS + ?i^c/2]/J, 



(38) 



where the cyclotron frequency is Wc = eB/m*, and Eqs 
and m* are the B = values of the dressed hole ground- 
state energy and effective mass, respectively. This predic- 
tion is in excellent agreement with the results where the 
Aharonov-Bohm interference is turned off. In contrast, 
the full results show additional quadratic i?-dependence. 
This is not so surprising, since for small magnetic fields, 
the effective hoppings depend quadratically on B through 
the Aharonov-Bohm interference terms, and one would 
expect that dependence to be mirrored here through the 
effective mass. 

We do not perform a more quantitative analysis 
because, as already mentioned, the A^ = 3 analytical 



case is not fully converged at this value oi t/J = 3; for 
example, the effective mass is m*/m « 54, whereas the 
converged result is m* /m sa 40 - this is quite a sizable 
difference. Going to a smaller t/J is not easy either, 
because the bandwidth narrows (or, equally, the effective 
mass increases) and it becomes more and more difficult 
to separate various features. 

V. SUMMARY 

To conclude, we have investigated the motion of a hole 
in a 2D square Ising AFM. We showed that summation 
of the contribution of all Trugman loops up to a given 
length N can be carried out numerically, and convergence 
is reached for A^ = 7 if i/ J < 3. Qualitatively correct 
and quantitatively quite reasonably accurate results are 
obtained from a simple analytical approximation for N = 
3. We find that the effective mass of the hole can be fairly 
low, of around 30. . . 40 m, if t/J ~ 3. . . 5. 

If a magnetic field is turned on, Aharonov-Bohm in- 
terference between clockwise and counterclockwise Trug- 
man loops leads to dependence of the effective hoppings 
on the magnetic field over and above the usual Peierls 
phases. Their most spectacular manifestation is in the 
"collapse" of the Hofstadter butterfly structure at fields 
where the Aharonov-Bohm interference is destructive. Of 
course, creating large enough magnetic fields to see the 
Hofstadter butterfiy band structure in a crystal is still 
a challenge (moreover, one would need a large J in or- 
der to prevent a transition to ferromagnet order, at such 
large fields). Nevertheless, this is an interesting effect 
that might be mimicked in some other type of systems. 

For small magnetic fields, where the magnetic length 
is large as compared to the unit cell. Landau levels form. 
Here, we find that the effect of the Aharonov-Bohm 
phases is to bring additional dependence on B through 
m*{B), in the cyclotron frequency. This is a small ef- 
fect, but one that might be easier to see through various 
quantum oscillations-type measurements. 
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